Triplet proximity effect and odd-frequency pairing in graphene 
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We study the interplay between proximity-induced superconductivity and ferromagnetism in graphene by 
self-consistently solving the Bogoliubov-de Gennes equations on the honeycomb lattice. We find that a strong 
triplet proximity effect is generated in graphene, leading to odd-frequency pairing correlations. These odd- 
frequency correlations are clearly manifested in the local density of states of the graphene sheet, which can be 
probed via STM-measurements. Motivated by recent experiments on S|N|S graphene Josephson junctions, we 
also study the spectrum of Andreev-bound states formed in the normal region due to the proximity effect. Our 
results may be useful for interpreting spectroscopic data and can also serve as a guideline for future experiments. 
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Introduction. Graphene^ constitutes a new exciting setting 
for studying the interplay between different types of long- 
range orders, such as ferromagnetism and superconductivity. 
Although the intrinsic appearance of both of these phenom- 
ena only occurs under special circumstances in graphene (see 
e.g. Refs. |2]Q, they can always be induced via proximity 
to host materials with the desired properties. The study of 
how the peculiar electronic properties 5 of graphene interact 
with superconducting correlations has recently attracted much 
attention both theoretically 6 8 and experimentallyJSE] Such 
studies are done by depositing two superconducting leads on 
graphene in order to create a graphene S|N|S Josephson junc- 
tion. By also exposing the N region to a ferromagnetic host, 
a hybrid S|F|S junction is constructed, which then will of- 
fer an excellent platform in which to study the interplay be- 
tween ferromagnetism and superconductivity. Previous work 
on such hybrid structures have reported on interesting effects 
by study ing it s transport properties via a scattering matrix 
approach! 12 ' 13 > However, this formalism does not include the 
full extent of the superconducting proximity effect, as it does 
not self-consistently solve for the superconducting order pa- 
rameter inside the junction. A self-consistent solution, on 
the other hand, will explicitly include the Cooper pair deple- 
tion in S, and the corresponding leakage into N, near the in- 
terfaces. In addition, the scattering matrix approach cannot 
deduce the symmetry of the induced correlations in the non- 
superconducting region, nor calculate their manifestation in 
the local density of states (DOS). In par ticula r, it remains be 
clarified how odd-frequency correlations 14 15 adapt to the un- 
usual electronic environment of graphene and what their sig- 
nature is in experimentally accessible quantities. The issue of 
odd-frequency pairing has recently generated much activity in 
the field of conventional F|S junctions,^ but has not yet been 
addressed in the context of graphene. 

In this paper, we present a self-consistent lattice-study of 
the interplay between ferromagnetism and superconductivity 
in graphene, with focus on the behavior of the local DOS to 
probe this interaction. Specifically, we investigate an S|F|S 
junction and the concomitant manifestation of odd-frequency 
pairing. To make contact with the current experimental status, 
we also consider an S|N|S junction and identify the appear- 
ance of subgap Andreev-bound states in the DOS. Our results 



should provide a guideline for spectroscopic measurements in 
graphene when superconductivity and/or ferromagnetism are 
induced by means of the proximity effect. 

Theory. We start out with the following tight-binding lattice 
Hamiltonian for graphene 

n = 1 E Vlss* + 9lh«) E iM*(flfu, + gin*) 
- E ^(/it/it/i/u + gl^glgii), (D 

and largely employ the notation and methods of Ref.[l7] Here 
fia Cffio-) is the creation operator on the A (B) site of the hon- 
eycomb lattice, t ~ 2.5 eV is the nearest neighbor hopping 
parameter, denotes summation over nearest neighbors, 
and a is the spin index. Moreover, = fii + ah\ is the spin- 
dependent chemical potential. We will assume that the native 
chemical potential \l\ is a constant within each region of the 
junction (S or N/F). Experimentally, an overall chemical po- 
tential can be set in the whole sample by applying a back gate 
voltage. In addition, it is expected that some charge transfer 
takes place between graphene and the superconducting leads 
and, therefore, fi$ can sometimes be higher than fij^/p. hi 
is the site-dependent exchange field, which regulates the fer- 
romagnetic order induced by proximity to the ferromagnetic 
host material. Thus, is only non-zero in the region between 
the two superconducting leads. The last term in Eq. ([T]i mod- 
els the influence of the superconducting leads on graphene. 
The attractive on-site interaction U\ gives rise to s-wave su- 
perconductivity and this parameter is only non-zero in the S 
regions of the junctions. 

We perform a mean-field approximation on Eq. Q, with 
the superconducting order parameter defined by 

A ; = -DiK/n/tf) + <<7u<7it>]/2- (2) 

We have ignored spatial variations within one unit cell and 
used Aj = A g . We further consider a geometry with trans- 
lational invariance along the interfaces, i.e. orthogonal to the 
direction of the junction. For concreteness, we focus on a zig- 
zag interface and Fourier-transform the eigenvectors in the y- 
direction. Please note that for an ,s-wave symmetry, the spe- 
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FIG. 1: (Color online). Local DOS on a large energy scale (left panel) and near the gapped region (right panels) in an S|N|S graphene junction 
for a variety of lengths and FVMs. Large values of the DOS are indicated by a bright color, while small values are given by a dark color. We 
have fixed L$ = 45 sites and fls/t = 0.6, as well as set the superconducting phase difference to zero, (j> — 0. (a): /Ltjv = Ms ( no FVM), and 
Ljv = 8 sites (£ ~ 25 sites), (b): fix = fis and Lm = 24 sites, (c): /xjv/i = 0.2 and Ln — 8 sites, (d): /j,m jt = 0.2 and Ljv = 24 sites. As 
seen, in-gap bound states, below the gap edge, are formed and persist even in the long-junction regime. 



cific direction of the interface will not matter. After diagonal- 
izing the problem, we arrive at the tight-binding Bogoliubov- 
de Gennes (BdG) equations for graphene 



fHf(n,m) A(n,m) 



^\aU \ fr I W C = ^(MC- (3) 
^ \A^(n,m) -H±(n,m)J 

Here n is the lattice site index along the junction and 



k y = 2x1/ ( N y a), where I is an integer such that k y € 
] — tv /a, x/a], and a is the lattice constant. Moreover, ip£ = 
[<{ky)> Vn(ky), <( k y)i <Os)F is generated by two copies 
(one for each sublattice) of the standard canonical Bogoliubov 
transformation. We have also introduced the following matri- 
ces 



H a (n,m) 



-t\S„ 



-(fin + ah n )S 

+ 26 n -i >m cos(k y a/2)] 



+ 2S n+lim cos(k y a/2)} 
-(fi n + ah n )5 

n id 



, A(n,p) 



A„(5„ 
Q 



(4) 



For a self-consistent solution of the above equations, we first 
guess an initial A n , then find the corresponding eigenvalues 
and eigenvectors to Eq. ([3]), followed by a recalculation of A„ 
using the self-consistency criteria in Eq. (|2]). This process is 
iterated until A„ no longer changes between subsequent iter- 
ations. We note in passing that for the case of a zero exchange 
field, Eq. Q is particle-hole symmetric and it is then enough 
to solve for only the negative eigenvalues. However, this is 
not the case for a non-zero h, and we are therefore forced to 
calculate all the eigenvalues. 

We are particularly interested in investigating the appear- 
ance of spin-triplet correlations in the system when the ex- 
change field hi is non-zero. These triplet correlations are 
necessarily of an unusual nature, since the spatial symme- 
try of the superconducting order parameter is isotropic (s- 
wave). According to the Pauli-principle, it is possible to 
have superconducting correlations which are both isotropic 
and spin-triplet simultaneously as long as these have an odd- 
frequency symmetry. The frequency dependence is obtained 



by Fourier- transforming the relative time-coordinate (t — t'), 
and it thus follows that such an odd triplet amplitude must 
be antisymmetric with respect to (t — r'). In effect, this 
amounts to a strong retardation effect since the correlator van- 
ishes at equal times. The study of proximity-induced odd- 
frequency pairing has recently generated much interest in con- 
ventional metallic S|F systems, but has not yet been explored 
in graphene. To investigate the presence of odd-frequency 
pairing, we introduceC^ 



(/it(r)/i4(0) + /H(r)/ it (0)). 



(5) 



As with the order parameter in Eq. d2}, we define the effec- 
tive odd-frequency correlator T\ as the average between the 
expectation values on the two sublattices described by the /- 
and g-fermion operators. We will also be concerned with the 
local DOS Afi(E) which in the low-temperature limit is ob- 
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tained via the charge density p as follows: 



i = / Mi{E)dE = J2(fUi° + algur) ■ (6) 

J — oo „ 



Results and Discussion. In what follows, we present a nu- 
merical and self-consistent solution of the above equations. 
Let us first consider the S \N\S case, shown in Fig. [I] To model 
experimentally relevant scenarios, we consider both short and 
long junctions in addition to the presence or absence of a 
Fermi-vector mismatch (FVM) at the interface. It may be in- 
structive to start with a reminder of the analytical result for 
the Andreev bound-state energy inside a junction, obtainable 
in the ultrashort-junction regime where L <C £, with £ be- 
ing the superconducting coherence length. This relation reads 
E = A [l - Dsin 2 (</)/2)] 1 / 2 , where D denotes the inter- 
face transparency and </> the superconducting phase difference 
between the two leads. For <j> = 0, it is seen that the bound- 
state lies right at the gap edge, independent of the the inter- 
face transparency. As a consistency-check, we have verified 
that we also obtain this result numerically when using a non- 
selfconsistent step-function profile of the superconducting or- 
der parameter. Let us now turn to the self-consistent treatment 
in Fig.[T] where we consider four scenarios for a junction. In 
(a), we have no FVM (i.e. fig = Miv) an d a short junction, 
Ljy = 8 sites to be compared with £ ~ 25 sites. The su- 
perconducting regions are chosen to be large, L$ = 45 sites, 
so that they act as superconducting reservoirs. As seen, the 
self-consistent solution for the order parameter only slightly 
shifts the bound-states inside the gap. This can be understood 
as a consequence of the proximity effect suppression of the 
order parameter near the interface. In (b), we increase the 
length of the normal region to Lm = 24 sites. In this case the 
bound-states now reside well within the superconducting gap, 
in contrast to the analytical prediction for ultrashort junctions. 
Note that the large value of ensures that these states do 
not pertain to some surface-effect, but that they penetrate into 
the entire N region. When turning on a FVM (/is > //jy) in 
(c) and (d), it is seen that the magnitude of the in-gap bound- 
states in the N region is strongly reduced, due to the reduced 
normal-state DOS in the N region, although the general en- 
ergy dependence is seen to be similar as in (a) and (b). 

We now turn to the case where the region between the su- 
perconducting leads is ferromagnetic via a proximity to a fer- 
romagnetic host material. The ground-state of an S|F|S junc- 
tion can occur for a superconducting phase difference cf> of 
either or n, so it becomes necessary to consider the free 
energy of the junction to correctly identify the ground-state. 
In Fig. [2] we consider a short junction with Lp = 8 sites 
and no FVM. Our results remain qualitatively the same also 
when including a moderate FVM (e.g. hf/ps = 0.8). How- 
ever, for a sufficiently low chemical potential in the F region, 
the DOS becomes too small to sustain any appreciable in-gap 
electron density, as was also seen in Fig. [T] In the experi- 
mentally relevant scenario, it is reasonable to expect that /i p 
is doped away from the Dirac point to support the presence 
of ferromagnetism, which is exactly the case considered here. 
In Fig. |2ja), we show the spatial self-consistent profile of the 
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FIG. 2: (Color online). Local DOS in an SjFjS graphene junction 
with Ls = 50 sites, Lf = 8 sites, h/t = 0.05, and fis/t = 
fip/t = 0.6 (no FVM). (a): Spatial profile for the superconduct- 
ing order parameter. The position of the interfaces are marked with 
black vertical lines. The inset shows the imaginary part of the in- 
duced odd-frequency correlations, which from bottom to top corre- 
spond to times tr = 1, 2, 3, 4, 5. (b): Local DOS on a large energy 
scale (left panel) and near the gapped region (right panels). Large 
values of the DOS are indicated by a bright color, while small values 
are given by a dark color. As seen, odd-frequency correlations give 
rise to a zero-energy peak in the DOS. 



superconducting order parameter and the generation of odd- 
frequency correlations in the F region (inset). The magnitude 
of the corresponding anomalous Green's function peaks in the 
middle of the ferromagnetic region, and penetrates a short dis- 
tance into the S regions. We have here plotted the imaginary 
part of T l since it couples directly to the DOS. 19 The experi- 
mental manifestation of such odd-frequency correlations have 
previously been discussed in conventional metals, in which 
case one expects an enhancement of the DOS at the Fermi 
levelfH0 i n pjg j2jt>) we demonstrate that the same signature 
applies for graphene: a zero-energy peak emerges and it is 
also flanked by additional in-gap bound-states. The ground- 
state for the parameters used in Fig. |2]was numerically found 
to be the 0-phase. 

It is well-known that the superconducting phase-difference 
(\> can be tuned actively via an external flux or current flowing 
through the system. We investigate in Fig. [3] how the local 
DOS in a S|F|S junction changes when the length of the F re- 
gion is increased (Lp = 24 sites) and compare specifically 
the 0- and 7r-phases. As seen in Fig. |3ja), the odd-frequency 
pairing correlations are clearly manifested for the 0-phase as a 
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FIG. 3: (Color online). Local DOS in an SjFjS graphene junction 
with L s = 50 sites, L F = 24 sites, fj, s /t = 0.6, and h/t = 0.06. 
(a): 0-phase, fx F /t = 0.6 (no FVM). (b): vr-phase, fi F /t = 0.6 
(noFVM). (c): 0-phase, fi F /t = 0.2. (d): 7r-phase, fi F /t = 0.2. 
Large values of the DOS are indicated by a bright color, while small 
values are given by a dark color. As seen, the odd-frequency pairing 
amplitude is most pronounced in the 0-phase, both with and without 
a FVM. 



strong zero-energy peak when there is no FVM. Going to the 
7t -phase in (b), it is seen that this peak is suppressed whereas 
the non-zero bound-states inside the gap are instead more pro- 
nounced. Upon introducing a strong FVM in (c) and (d), we 
see again how the magnitude of the proximity effect in the 
F region is severely suppressed, although some signs of the 



odd-frequency correlations are still visible in (c) and (d). 

The presence of significant in-gap DOS in the whole F re- 
gion, even for long junctions (Lp = 24 sites), rules out the 
possibility of these being caused by surface-states at the S|F 
interfaces. We have performed numerical calculations for sev- 
eral sets of parameters to investigate the robustness of the 
odd-frequency peak, and find that it in general competes with 
the singlet correlations which instead induce a standard mini- 
gap in the electronic spectrum inside the junction. In spite of 
this coexistence, our results above demonstrate that the odd- 
frequency amplitude can be read out from spectroscopic in- 
formation in a feasible parameter regime. 

Summary. In summary, we have investigated in a self- 
consistent manner the proximity effect and its implications for 
the local DOS in both magnetic and non-magnetic graphene 
Josephson junctions. We have considered several experimen- 
tally relevant ranges of doping levels and junction lengths. It 
is found that a considerable triplet proximity effect can be in- 
duced in an S|F|S graphene junction, giving rise to so-called 
odd-frequency correlations. These are manifested clearly as 
zero-energy peaks in the DOS, which may be probed by STM- 
measurements. We have also identified the appearance of 
Andreev-bound states in a S\N\S graphene Josephson junc- 
tion. They appear as in-gap resonances in the DOS, in con- 
trast to non-self-consistent results positioning them at the gap 
edge. Our results should be helpful for the interpretation of 
spectroscopic data and will hopefully serve as a guideline for 
future experimental activity. 
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